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Abstract 

The paper is dealing with semi-classical asymptotics of a characteristic 
function for a stochastic process. The main technical tool is provided by 
the stationary phase method. The extremal range for a stochastic process 
is defined by limit values of the complex logarithm of the characteristic 
function. The paper also outlines a numerical method for calculating 
stochastic extrema. 

1 INTRODUCTION 

The extremum for a stochastic process admits transparent numerical presen- 
tation in terms of limit set of its characteristic function which is treated as a 
high frequency integral [7], [10], [II]j [12] • The proposed concept of stochas- 
tic extremum is compatible with other known methods of assessing extrema of 
stochastic functions (see, e.g., [5], [3], 0], [5], [S]). However, the ideology 
of high frequency integrals, though different, is close to the simulated annealing 
technique [5] , [3] ■ In our approach the role of the parameter like the inverse of 
"temperature" (from the annealing process ) is played by the frequency and in 
order to calculate the extremum we increase the frequency. Method of high fre- 
quency integrals is able to calculate extrema of stochastic processes with several 
variables, and therefore can be applied to analyze images and three-dimensional 
data samples. When it is possible to apply the method of assessing extrema with 
principal component functions (in the sense of Karhunen-Loeve representation) 
[6], then our approach works as well and practically leads to the same results. 
However, it does not rely on any type of Karhunen-Loeve representations and 
in this sense our concept of stochastic extremum is of equal or more general 
nature than the Karhunen-Loeve representation itself. Moreover, the method 
of treating stochastic extrema as stationary phases of its characteristic function 
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leads us to a transparent numerical procedure that allows efficient estimation of 
the stochastic extremum and evaluation of its statistical significance. 



2 STATIONARY PHASE 

Our goal is to introduce the concept of extremum for a stochastic process. In 
order to do that we employ the high-frequency integrals: 

/>oo 

I{k,uj) = / ip{t,uj)e'^-f^'^'^'>dt, 



where uj € is a, parameter; both / and ip are real functions that are infinitely 
many times differentiable with respect to t. Moreover, ip{t, lu) has a finite time 
support for any fixed lo G fl, suppt{ip) is a subset of a closed interval from 
IR. Throughout the paper IR denotes the set of real numbers. For reader conve- 
nience, we recall the basic properties of the high-frequency integrals (for further 
reading on this subject see, e.g. [10], [H], [12]). 

If suppt{ip) C [a, b] and 

— f(t,uj)^OVte\a,b] and V ui e n 
dt 

then 

I{k,uj)^ (p{t,uj)e'''-f^'^'^'>dt= (p{t,uj)e'''-f^*^'^'>dt 



and integrating by parts n times yields 



J a 

where the linear operator L is defined as 

and ft denotes the derivative of f(t,uj) with respect to time, 

/t(i,c.) = ^/(t,c.). 

As one can see, if the phase f{t,Lu) does not have critical points in suppt{ip) 
then 

l(k,uj)^0(^) ynelti and ^uj e n, 

k" 

where IW denotes the set of natural numbers. This fact is often represented as 

l{k,u;) = 0{-^) Vioen. 

The critical (or stationary) points of the phase f{t,uj) make the main contribu- 
tion into the high-frequency integral I{k,Lu) as k —^ oo. 
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Definition 1 A point {t*,u}*) gK x Cl is called a stationary phase (point) if 

The set of all stationary phase points for / is addressed as St{f) C IR x Q. 
Now let us turn our attention to calculating the contribution of a stationary 

phase point {t*,uj*) € St{f) into the high-frequency integral I{k,uj). A station- 
ary phase point {t*,cu*) G St{f ) is said to have order m G W if m is the first 
natural number for which 

The set of such stationary phase points is denoted by St'^{f). The Taylor ex- 
pansion near t* is 

fit, CO*) - fit*, uj*) = uj*)it - t*r + out - t*r+') as t ^ t*. 

Consider the change of coordinates 

xit,u;*) = isignif^^\t\oj*)) ■ ifit,co*) - fit\u;*)))^, 
where signiff"^\t* ,lo*)) denotes the sign of 

i±rfit*,uj*). 

Since 

±xit,co*) = \±ilrfit*,uj*)\^+Oit-t*) 
at ml at 

the change of coordinates is not degenerate on some interval Q^, 

t*-e<t<t* + e 

Let us take infinitely many times differentiable function h with supptih) C Qe 
(the set of such functions is denoted as Co^iQe))- Assume also that hit, w*) = 1 
in a neighborhood of t*. Then 



Ut 



ifit, uj*)hit, uj*y^-f^^^'^*Ut+ / (^(i, oj*)il-hit, uj*))e'''-^^*''^* 

-OO J — oo 

and in order to find the contribution of the stationary phase (i*, cu*) we need to 
calculate asymptotics for 

/oo 
>fiit,LO*)hit,u*)e'''-f'^'^'^*Ut. 
-OO 

After making the change of coordinates x = x{t^uj^) in the integral 

/oo 
<pit, uj*)hit, w*)e^'=(/(*''^*)-^(**''^*))df 
-OO 
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we have 

Iik,Lu*) = e''=/(**^-*) . r ^{x, Lu*) ""^^ e-'g"(^^' dx, 
where a:* denotes ^x{t,uj*). The integral 

/>CXD 



can be calculated by reducing it to the linear combination of the integrals like 



and then evaluating the latter with the help of the integral along the curve on a 
complex plane [12j . The curve consists out of the segment of x-axis < a; < p, 
the arc of a circle 

p-e±- (0<r<^) 
and the segment of the straight line 

r-e=^'^ (/9>r>0). 

Taking p — > oo yields that 

/ Q±ikx (^2; = cos( )^--C„i for odd m 



and 



e'^^''^ dx = - — J — -Cm for even m 
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oo 



where 

Cm 2 

Taking into account that h{t'*,uj*) — 1, for even m we have 

I{k, uj*) = ^{t*, Lj*) ■ C„, ■ ( ^ ) " • e"^3"(^*'"' . e''^-^^** + 

^ • \ ft 

and the latter integral has the asymptotic 

0{ — 2") as fc ^ cxD 
k ^ 

as long as there are no other stationary phase points present. If m is odd then 
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as A; — »■ oo. 

In conclusion to this section we formulate the basic results of stationary 
phase method in the form of the following formal statement. 

Theorem 2.1 Let (p{t,u!), f{t,uj) be real functions that are infinitely many 
times differentiable with respect to t for any w G O. Moreover, for any w G O 
one can find an interval [a{u)), b{co)] C IR such that 

ip{t,oj)eC^{[a{w),b{w)]). 

Then the following statements hold. 

i. If [a{w),b{w)] n St{f) = then 

I{k,u}) =0{-^) as k^oo 

a. If [a{Lv), b{oj)] n St{f) = {tj} then I{k,uj) has the asym,ptotic 

even mj \ l^ ' \Jt V-ji^)] 

k"-' / Oddm, V ^■\ft 



as k —> oo. 



3 STOCHASTIC EXTREMUM 

Consider a real valued stochastic process ^{t) defined in the probability space 
(ri, A{U), P), where A{n) is a cr-algebra of subsets from and P is a measure of 
probability. Throughout the paper we assume that $, {t) takes only non-negative 
real values. In the context of this paper it is tacitly assumed that 

^{t,uj) = g{t,u!), 

where a; is a stochastic variable which docs not depend on time and g{t, x) is a 
smooth non-negative function of its arguments. 

In order to examine whether ^{t) has extrema on interval [a, b] C IR we take 
a function ips{t) e C^{[a -s,b + s]) such that s > 0, ips{t) > V t G IR and 

<^c(t) = 1 Vt G [a,b]. 
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Then we analyze the asymptotics of the high-frequency integral 

Jq J -oo 

The formal description of the situation when ^{t) does not have stochastic ex- 
trema on [a, h] sounds as follows. 

Definition 2 Let ln(z) denote a fixed branch of the complex logarithm. Then 
a real-valued stochastic process ^ (t) does not have stochastic extrema on [a, b] if 
one can find a positive real number e such that 

lim Re{\ ln( / / (p,{t)e'''^^*''^UtdP{uj))} = 0. 
fe^oo tk Ja J_oo 

The logical negation of this statement describes the intervals where stochastic 
extrema for ^(t) occur. In other words, the interval [a, b] contains stochastic 
extrema if V £ > 

lim Re{\- ln( / / ip,{t)e'''^^*''^UtdP{uj))} 0. 
fc^oo tk JnJ-oo 

Sometimes it is possible to estimate extremal values for ^ (t) on [a, b] with the 
help of the following 

SmaXe{[a,b]) = ]h^k^ooRe\^H[ [ (p^{t)e'''^^^''^^dtdP{uj))\ (1) 



and 



Smine{[a,b]) = ljmk^^Rei^^\n{J J >fie{t)e'''^^''^^ dtdPiiv))^, (2) 

where lim and lim denote upper and lower limits, respectively. 

For majority of applications SmaXs ( [a, b] ) , Sm,in^ {[o-ib]) deliver boundaries 
for extremal values of a real stochastic process ^{t). One can justify that when 
the stochastic process ^{t) has additional properties Uke, for example, is a 
smooth manifold and ^{t,uj) is a smooth function of t and w. Let ^ denote 
the probability density function for P{u)). Then the following statement takes 
place. 

Theorem 3.1 Let be a smooth manifold. Assume also that density ^ ex- 
ists; ^ e C°°{fl), ^{t,uj) e C°°(IR,r2) and there is only one stationary phase 
point {t*, u!*) e St'^{^) n{[a - e,b-\- s] x O) with ^^(i*, w*) = and the second 
derivative has non-zero determinant 

J{t\oj*)=det(-^{t\u;*)\ j^O (3) 
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at the stationary phase point {t* , lu*). If ■^P{lu*) > then the following is true 



ik 



Lp,{t)e'^^^'^'^UtdP{Lo)j^ ^ i{t\L0*)+O{^) as k 



Proof. 

Let n denote the dimension of CI. Then applying the stationary phase method 

m M to 

/• /'OC 

if,{t)e'''^^'^'^^dtdP{Lu) 



yields 

/• /"OO 

ipe{t)e'''^(*UtdP{uj) 



1+1 



^^ >^ — OO 

duJ ^ ' k' ,/\J{t*,LO*)\ 

where 

denotes the difference between the number of positive and the number of neg- 
ative eigenvalues of the corresponding quadratic form. Taking the complex 
logarithm of 

A. pf^^Af^^T^ . ^ . e^^9»(iSlr(**."*))*f . ''^'') . fl + 0(-)) 

duo ^ k' ^\J{t\uj*)\ ^k" 

we obtain 

^,{t)e'^^^'Ut) = iki{t\oj*) _ln(fc)+ln(4-PK)) + ^^-ln(27r)- 

-oo 2 2 

\ In (|J(^^ .'^)|) + ■ ^^9n{^{t\ u*)) + 0(1) 

as k ~* oo. Dividing by ik and taking real part complete the proof. Q.E.D. 

Formulas H]), ([2]) together with Theorem 13.11 give us a recipe for assessing 
extrema of a real stochastic process ^{t) on an arbitrary time interval [a,b]. 
In order to do that one needs to analyze the limit properties of the following 
integral for big values of k : 

ikj^^ z 

where the complex curve 7^ is defined as 

7fc = ( / / ^,(Oe'^«(*'")dtdP(w); < A < fc 
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Notice that 7^ can be interpreted as the characteristic function of ^(i, uj), where 
V5e(t) is chosen so that 

/oo 
(p^{t)dt = 1. 
-00 

In conclusion we notice that the condition ^ can be relaxed with the help of 
Theorem 12. II fscc also [T], [S]). In some applications (e.g., when estimating tails 
of distributions ) it is beneficial that the truncation function ipi^ (t) also depends 
on UJ. 
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